XA0-A135  371 


IDENTIFICATION  OF  A  CLASS  OF  FILTERED  POISSON  PROCESSES  1 // 
(U)  NORTH  CAROLINA  UNIV  AT  CHAPEL  HILL  DEPT  OF 
STATISTICS  D  DE  BRUCQ  ET  AL .  1981  N00014-75-C-0491 

F/G  12/1 


UNCLASSIFIED 


NL 


one  file  copy 


TITLE  :  IDENTIFICATION  OF  A  CLASS  OF  FILTERED  POISSON  PROCESSES 


Authors  :  DE  BRUCQ  Denis*  -  GUALTIEROTTI  Antonio** 


Summary  :  The  formulation  of  a  model  is  an  essential  step  in 
the  experimental  process  of  understanding  a  phenomenon.  In 
this  paper  a  new  class  of  filtered  Poisson  processes  is  intro¬ 
duced  :  the  amplitude  has  a  law  which  is  spherically  invariant 
and  the  filter  is  real,  linear  and  causal.  It  is  shown  how 
such  a  model  can  be  identified  from  experimental  data. 


* 

Professor  University  of  Rouen  -  BP  67 
76130  MONT  SAINT  AIGNAN  -  France 


** 

Professor  IDHEAP,  bdtiment  des  Facultds  des  Sciences 
Humaines  I  -  University  de  Lausanne 


CH-1015  -  LAUSANNE 

Research  supported  by  ONR  contracts 

N00014-75-C-0491 

N00014-81-K-0373 


DT1C 

ELECTE 
DEC  6  1903 

D 


distribution  statement  a 

Approved  I  or  public  ftlfciMi 
_ Distribution  Unlimited 


83  12 


078 


3 CCU HlTY  CLARIFICATION  OF  THIS  PAGE  (Whom  Of 


REPORT  DOCUMENTATION  PAGE 


*.  TITLE  (m*  MMIIt) 

Identification  of  a  Qass  of  Filtered  Poisson 
Processes 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


a.  recipient's  catalog  number 


3-  Type  Of  REPORT  *  PERIOD  COVERED 

Technical 


«•  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR^) 


Denis  De  Brucq  5  Antonio  Gualtierotti 


3.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina 


11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


Statistics  §  Probability  Program 
Office  of  Naval  Research 

Arlington,  VA  22217  _  _  _ 


MONITORING  AGENCY  NAME  A  AOORESSfi/  different  trom  Controlling  Of  tic*) 


••  CONTRACT  OR  GRANT  NUMtEROI 

N00014-75-C-0491 

N00014-81-K-0373 


10.  PROCRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

NR  042-269 
SRO  105 


I*.  REPORT  DATE 


<3.  NUMBER  OP  PACES 

39 


15.  SECURITY  CLASS,  (of  thlm  report) 

UNCLASSIFIED 


tSa.  OECLASSIHCATIQN/ DOWNGRADING 
SCHEDULE 


1*.  DISTRIBUTION  STATEMENT  (ol  thi*  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  tn  Utock  20,  It  <u/t*r*nt  trom  Report) 


13.  EEV  BOROS  (CoMlnvm  on  reyetee  mlde  It  nocoeemry  "nd  Identify  by  block  number) 


Filtered  Poisson  processes;  spherically  invariant  processes;  wave  propagation 


20.  ABSTRACT  (Continue  on  roeeree  *tdo  It  nocoommrf  ond  Identify  by  block  number) 


The  fomulation  of  a  model  is  an  essential  step  in  the  experimental 
process  of  understanding  a  phenomenon.  In  this  paper  a  new  class  of  filtered 
Poisson  processes  is  introduced:  the  amplitude  has  a  law  which  is  spheric¬ 
ally  invariant  and  the  filter  is  real,  linear  and  causal.  It  is  shown  how 
such  a  model  can  be  identified  from  experimental  data. 


00  \  jmTt*  1473  COITION  OF  1  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED  ~  ' 

SECURITY  CLASSIFICATION  OF  THIS  RAOS 


fFVitn  Brno  ImmnQ 


i 


INTRODUCTION 


The  process  introduced  in  this  paper  serves  as  a  useful  model 
for  the  study  of  random  wave  propagation  problems  and  certain 
areas  of  electronics  and  biochemistry.  An  initial  phenomenon, 
modeled  as  a  Poisson  point  process  whose  effect  is  described 
by  a  spherically  invariant  random  variable  C  is  observed 
through  a  causal  linear  filter  whose  impulse  response  is  G. 

The  observed  process  has  the  form 

ao 

z(w,  t)  =  l  c.(cd)  g ( t  -  -r .  (co> ) ,  -  -  <  t  <  +  • 
j-»  J  J 

where  Tj  is  the  time  of  the  j  the  jump  of  the  original  Poisson 
process . 

Spherically  invariant  laws  are  mixtures  of  normal  ones 
^ a '.■<?  y  in  i-he  st.tidv  of  wave  pronacat 1 on 

problems 

The  first  section  of  the  paper  contains  a  characterization  of 
such  laws  based  on  Choquet ' s  representation  theorem  [4]  and 
some  examples  (the  exponential  and  the  student  probability 
laws) . 

In  the  second  part  of  the  paper  one  finds  the  characte¬ 
ristic  function  of  the  process  Z  defined  above,  which  is 
strongly  stationary.  This  is  used  in  the  third  part  to  show 
that  the  odd  moments  of  Z  vanish  and  to  compute  the  moments  of 
order  two  and  four. 

The  identification  of  the  process  Z  requires  that  the 
impulse  response  G  be  estimated  knowing  the  moments  of  Z.  This 
can  be  done  if  the  filter  is  assumed  to  have  minimal  phase,  a 
concept  frequently  used  in  automation  (cf  [6] ) . 


y 


The  fourth  part  of  the  paper  presents  an  algorithm  for 
approximating  G  based  on  the  Fourier  transform. 

The  final  and  fifth  part  of  the  paper  is  devoted  to  the 
estimation  of  the  parameter  A  of  the  original  Poisson  process 
and  of  the  mixing  law  u  which  yields  the  spherically  invariant 
one  used  in  the  model .  Here  again  estimation  is  by  the  method 
of  moments  :  the  estimator  of  A  is  obtained  on  the  solution  of 
a  linear  equation  and  to  estimate  u  one  assumes  it  is  a  convex 
combination  of  point  masses  and  adjusts  the  classical  solution 
of  the  moment  problem  to  the  case  of  a  support  contained  in 
the  positive  half-line. 
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We  wish  thank  Professor  C.  R.  BAKER  of  the  University 
of  North  Carolina  at  Chapel  Hill  (Department  of  Statistics) 
for  many  helpful  discussions. 


I  -  CHARACTERIZATION  AND  EXAMPLES  OF  SPHERICALLY  INVARIANT 


RANDOM  VARIABLES 


The  problem  considered  here  is  that  of  determining  a  class 
of  spherically  invariant  distributions  (see  [5])  on  the  real 
line.  The  extreme  points  of  this  set  are  the  Gaussian  laws  Pg 
with  density  : 


(1) 


x2 
2a  2 


where  dx  is  the  Lebesgue  measure  and  one  wants  to  characterize 

the  closed  convex  set  generated  by  the  family 

{Po,  o  >  0}.  In  terms  of  characteristic  functions,  one  wants 

to  solve  the  equation  : 

a2  u2 

°°  —  -  ■  -  -  +ao 

(2)  $ (u)  =  j  e  2  u (da)  =  /  eiux  P(dx) 


where  <j>  is  the  characteristic  function  of  a  spherically  inva¬ 
riant  distribution  P  for  u  if  t  is  given  and  for  $,  if  ii  is. 


I . I  -  Lemma  : 

If  P  is  spherically  invariant*  P  is  symmetric  with  respect 
to  the  origin. 


P*oo  4  :  If  P  is  spherically  invariant,  then  by  definition  or 
reference 


4*oo 


+(u)  =  /  eiux  P (dx)  =  /  e  *  u (da) 

—  oo  O 


Since  <Mu)  *  ♦(-  u),  P  is  invariant  under  the  transfor¬ 
mation  x  -  x.  ■ 


*  Some  author  takes  a  larger  class  than  the  one  considered  here 
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We  know  that  <M0)  *  1  (then  u(iR+)  =  1)  and  that  t  is 
continuous  in  a  neighborhood  of  zero. 

1.2  -  Lemwa  : 

If  P  is  spherically  invariant,  there  exists  in 
C*(]0,  •  [)  (with  derivatives  of  any  order  in  ]0,  <■[)  such 
that  for  every  real  u 

(3)  *(u)  =  Mu*) 

P\oo&  :  Define  41  by  the  relation 


Cv)  ^  /  e 


du  (a) 


Then  for  every  positive  v  and  natural  integar  n 


/  -^ne‘^V 


du (a)  <  » 


so  that  4*  may  be  differentiated  arbitrarily  often  in  order 


to  obtain  : 


«*"  l. 


n  -  v 


du (a)  ■ 


The  transformation  defined  by  T  :  a+  b  ^  is  a 
bijection  of  fil  +  .  Let  v  »  u  o  T~l,  the  image  v  of  u  by  T. 

Then  *(v)  *  /  e”b  v  dv(b) 

«U 

1.3  -  Theorem  :  (  [4}  p.  237)  : 

There  exists  a  positive  measure  v  onl(t  such  that 

V  v  >  0  4»(v)  ■  /  e“b  v  dv(b) 


/ 


-  7  - 


if  and  only  if 

\\>  is  real  valued  on  IK  ,  is  ^  and  1  n  ^  0  (-  1 ) n  — -  >  0 

dvn 

1.4  -  T heo^em  : 

P  is  spherically  invariant  on  OK.'fl)  if  and  only  if 

(4)  P  is  symmetric  with  respect  to  the  origin 

(5)  the  characteristic  function  <j>  of  P  belongs  to  fe"(fiV\{0}) 
and  the  function  ,  defined  by 

(6)  ^ut'K  <f>  (u)  -  i|/ (u2 )  is  such  that 

VnfN  Vv  >  0  (-  1)n  (v)  >  0 

dv 


PtiooJ  :  If  (4),  (5)  and  (6)  hold,  one  has,  by  1-3,  that 

(7)  Mv)  =  f  e-b  v  dv  (b) 

This  can  be  written  as 

_  a2  u2 

— 5 — 

(u)  =  4»  (u2 )  =  /  e  du(a)  which,  as  we  know  is 

IK* 

the  characteristic  function  of  a  spherically  invariant  proba¬ 
bility  . 

If  P  is  spherically  invariant,  (4) ,  (5)  and  (6)  follow 
from  lemmas  1  and  2 .  B 

1.5  -  Example  : 

The  double  exponential  j  e“^x^  is  spherically  invariant. 
Indeed  a  straightforward  calculation  yields  : 


/ 


-  8  - 


a‘u* 


a*  XJ 


iux  X  -X  x  , 
e  2  e  '  1  dx  = 


X2  +  uJ 


-  / 


:he  result  is  obtained  by  settina 

a2  X2 

g—  (a)xx2  a  e  2  or 


4>  (v)  = 


X2  +  V 


and 


X 2  a  da 


an 


dv 


<P  (v)  =  (-1)  (-2)  ...  (-  n) 


(X2  *  V) 


n  +  1 


1.6  -  Example.  : 

X (e“au2  -  1 ) 

<t>  (u)  =eA  X  >  0  a  >  0  is  the  characteristic 

function  of  a  spherically  invariant  function. 

a 2  i 

Define  T  by  a  >  ^  and  v  by  v  =  u  o  T 

Since  T  is  a  bijection  of 
a2  u2 

ao  —  -  ao 

<Hu)  =  J  e  2  u  (da)  =  /  e_boU  v(db) 
o  o 

Choose  for  v  a  Poisson  law  with  parameter  X.  Then 

4>(uJ  =  /  e  2  4(da)  =  eX  ~  ^ 

o 

1.7  -  Example  : 

Student's  law  with  n  degrees  of  freedom  is  spherically 
invariant. 

y 

Student's  law  is  the  law  of  the  ratio  ^  of  two  indepen¬ 
dent  random  variables  X  and  R ,  the  law  of  X  being  normal  with 
mean  0  and  variance  a1  and  that  of  R  being  a  with  n  degrees 
of  freedom. 


Thus  the  joint  density  f  of  X  and  R  is  given  by  : 
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1 


f(x,  r)  =  -  exp  - 

/2n  2a 2 


D.  _  •, 

r(|)  2  2 


r*  r  1 

exp  -  Z—  (^)  — 


2a 


a  o 


The  characteristic  function  <|>  of  —  can  be  computed  di- 

rp  K 


rectly.  By  definition  one  has 


.  x 

iu  o 

y  u  *  'X  4>  (u)  =  E  (e  R) 


iu  ^ 


-  /  /  e  r  -L.  exp  -  2iL_ 
/27t  2a2 


exp 


n 


r(f)  2  1 


-  1 


r*  .r«  n~^  dx 
,2  °  a 


2oJ 


Applying  Fubini's  theorem,  one  gets,  integrating  over  x. 


4>t(u)  *  / 


2r 


r(§>  2* 


-  1 


2  r 

(§) 


n-1 


2a2  dr 


To  obtain  the  usual  form,  of  the  characteristic  function 
set  a  ^  p  ;  this  transformation  is  a  homeomorphism  of 
]  ') ,  » £  so  that 


4>m<U)  * 


J 

^2 


1 

2a2 


r(|) 


da 

n+1 


-  -  1 


f  r 
I 


/ 


-  10 


II  -  THE  CHARACTERISTIC  FUNCTIONAL  OF  THE  FUTEREV  POISSON 
PROCESS 

Let  C  be  a  spherically  invariant  random  variable.  We  have 
seen  that  its  characteristic  function  has  the  representation  : 

a*  u2 

(1)  4>(u)  =  E(eiuC)  =  /  e  2  du(a)  where  u  is  a  proba- 

*♦ 

bility  measure  on  (it\+,  ft>  +  ) 

Let  r  j  be  the  time  of  the  j-th  jump  of  a  stationary 
Poisson  process  with  parameter  X . 


As  is  well  known  for  any  fixed  interval  [a,  b]  the 
number  N  of  events  t  .  in  [a,  b]  is  a  Poisson  random  variable 
with  parameter  (b-a)  so  that 

(2)  PIN  =  n)  =  e“X (b-a)  £x (b--a)j 


II  .  1  -  PfL0£eAtl£  i 

When  N  =  n  is  fixed,  the  times  ti,  tj,  ...»  tN_n  of 

intervals  [a,  b] ,  are  random  and  can  be  chosen  independent  and 

of  same  uniform  law  :  • 

b— a 

Finally,  a  jump  at  time  t,  of  normalized  sized,  has  an 
effect  described  by  G(t,  t) . 

If  we  suppose  the  stationarity  then 
G(t,  t)  =  G(t  -  t) 

If  we  suppose  the  causality  principle  then  * 


I 

l 


G(t-x) 


T  $  t 


G  is  called  the  Green  function  of  the  process . 


The  filtered  Poisson  process  is  then  defined  by  the 
relation  : 

(4)  V  tfft  Z(t)  =  l  C.(w)  G  ( t  -  t.) 

j<  £  3  J 

The  Poisson  process  and  the  amplitudes  Cj  ]  (-  &  are 
supposed  to  be  independent . 

In  what  follow,  the  following  convenient  notation  shall 
be  used  : 

(5)  ^TfiK  GT(s)  =  G  ( s) 

1  ^-j  being  the  indicator  of  the  interval  £0,  t3 

1 1. 2  -  P^opoa-c-t-cow  ••  If 

(6)  E ( | C | )  <  -  and 

oo 

(7)  /  |G(s) |  ds  <  »,  then 

O  oo 

Z(t)  is  integrable  and  E(|Z(t)|)  $  XE(|C|)  /  |G(s) |  ds 

o 

Piooj  :  Since  |2(t)|  $  \  |C^|  |G(t-r^)|»  one  has 


I 

-  12  - 


(8)  E  (  |  Z  (t)  )  |  x<  E  (  |  C  |  )  E(  l  |G(t~r.)|) 

j  **  3 

An  upper  bound  for  E(  \  |G(t-x.) |)  is  obtained  by  consi 

jf-4.  3 

dering  finite  time  effects  represented  by  the  expression  : 


(9)  E (  l  |G  (t  -  x  • )  | ) 

■*£  £  J 

GT<t  -  t^)  is  zero  except  for  Tj  in  [t  -  T,  t] ,  so  that  condi¬ 
tioning  on  the  number  of  Tj'  s  in  [t  -  T,  t]  one  has  using 
probabilities  (2)  and  property  (3)  : 

E(  l  |G  (t-T,)|)  =  l  e'XT  (  l  E(|G  (t-T  )  |/N  =  n)  ) 

jea  3  n=1  n*  j  =  1  1  3 


e 


-AT 


OD 


l 

n=1 


( AT)  n 
n ! 


t 


n  J 

t-T 


|  G (t-x ) 


dr 

T 


=  X  /  |G(s) |  ds,  so  that 

o 

T 

YTfli  E  (  l  |G  (t-T.)  |)  =  A  /  |  G  (s)  |  ds 

jfi  3  o 

Taking  the  limit  as  T  +  »,  one  has  the  needed  result  : 


(10)  E(  l  | G ( t-r . ) | )  =  A  J  | G (s) |  ds 
jfc*  3  o 


We  are  going  to  characterize  the  temporal  law  of  process  Z. 
To  have  a  practical  expression  for  the  characteristic  function  of 

(11)  Wn  *tx,  t,,  tmfik  Z^  (Z(ti),  Z(ta),  z(tm)) 

further  assumptions  are  needed,  as  follows  : 

OO  00 

(12)  J  |G(s)|*  ds  <  «  and  even  /  | G ( s) j 2  ds  =  1  and 

o  o 

00 

(13)  /  a*  du(a)  <  • 
o 

Since  E(|c|*)  =  J  du(a)  /  c*  — - —  exp  -  j  —  dc  *  /  a*  dn(a)  <  • 
*  *  a  1  a* 


(13)  means  that  C  has  finite  variance. 


s  2  2  » 

II . 3  -  Lemma  :  /  | 1  -  exp  -  Ga (a) |do  $  f  G* (a)  da 

o  o 

A  a  ^ 

PtooJ  :  Consider  f <£)  *  1  -  exp  -  i  -y  G2 (o)  then 
||  U)  =  G2  (a)  exp  -  l  G2  (a)  and 

i,  2  2  3 

/  G2 (a)  exp  -  A  G2 (a)  dX  *  1  —  exp  -  i  G2 (a) 

At  last,  we  take  i  =  1  then 
1  -  exp  -  G2 (o)  $  J  ^  Ga (a)  exp  -  X  G2  (a)  dX 

$  J  2i  g2  (a)  dX  =  G2  (a)  ■ 

O  * 

1 1 . 4  -  Theorem  : 

If  Z  is  the  derived  Poisson  process  : 

(4)  Ytfift  Zlt)  =  l  C.  G (t  -  x.) 

j*4  3  3 

if  (14)  and  (15)  are  satisfied  : 

|  |G(s) |2  ds  <  ®  and  J  a2  du(a)  <  »  then 

i\  + 

V  mfW  V(tx,  t, ,  ...»  tm>  t  iK  m  iu  =  (u2 ,  Uj ,  ...»  u^)  £ 
(14)  *(  u)  £  E(ei<u'Z  >)  = 

r  j  ®  * 

exp  -  X  /  ds  dix(a)  Il  -  exp  {-  ^-  (  J  G(tk+s))2} 
\k.il,  kal 

P*oo_£  :  Without  loss  of  generality,  one  may  assume  that 

ti  <  ta  $  ...^t  *  0  and  that  T  >  t  -  ti 

m  m 


I 
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-T  t 1  t2  t, 

Tl  T, 


Tj  T 


m-1 

—a  —  f —  x 


0=t 


m 


n-1 


By  definition  of  the  inner  product  and  of  Z ^ , 

(15)  <u,  Z^>  ft  ua  Z(ta)  +  ...  +  um  Z(tm) 

*  I  ci  l  uk  G(t  -  T.) 

jez  3  k-i  K  K  3 

as  done  previously  the  calculation  is  performed  first  for  finite 
time  effect  by  conditionning  on  the  number  of  jumps  : 

(16)  <u,  ZT>  »  l  C.  I  u.  G  (t.  -  x.) 

&  jei  3  k=1  K  T  K  3 

and  we  compute  the  characteristic  function  of  this  random  variable. 

T 

i<u,  Zr> 

(17)  ^T(u)  A  E(e  6  ) 

n 


=  e 


«  n  n  1  ci  I  uk  GT(tk'Ti) 

-XT  .  j  3-XT  -f  E(e  3  k-1  T  k  3  /H  .  „) 

nil  n!  j  =  i 


=  e 


•XT 


I 


V  <£  -V  V1* 


-vl" 


0  -  V  ^  Uk  GT(Vt)) 

_  +  X  /  J  drdu(a)  e  *  k  K  * 

=  e  e  -T  + 


=  exp  {-  X  J  J  drdu(a)  £l  -  exp{-  (1  uk  GT  (tk+T ) ) 


We  will  transform  the  exponent  of  the  exponential  ;  let  : 

&  rT  (  r  WVs))* 

(18)  *_(u)  S  log  ♦  (u)  -  -  X  /  /  dsdu(a)  1  -  e  *  k  K  T  K 

T  0  1^  L 


/ 


15  - 


and 


m 


(19)  G^(s)  =  I  uk  GT(tk  “  s) 

Then  by  lemma  I I. 3  : 

(20)  /  dsdu(a)  [l  -  exp  {-  y  Ga(s)  }]  $  {  dsdu(a)  y  G*  (s) 

*i  L  J  Hi  5 


m  m 


$  J  dsdn(a)  y  I  u*  1  GT(t,  +  s)a 
ik2+  2  k=1  *  k=1 


m 


m 


4:  /  dsdu(a)  y  £  G  ( t,  +  s )  *  \  u*  <  » 

FR^  2  k=1  *  k=1 


So  by  the  dominated  convergence 


T  00  r  2  m  1 

li***  f  f  ded.Li  * a '  •  1  —  pyp  ( —  — —  (  Y  n.  a  ( t.  .  +  s) )  *  1 ! 

TOO  L  *  k=1  K  1  K  J 

=  J  J  dsdu(a)  f"l  -  exp  {-%-(!  u.  G(t.  s) )  *  }1 
0  0  L  *  k  =  1  J 


So  we  obtain  : 


I  a2  I 

<t>  (u)  =  exp  -  xf  dsdu(a)  1  -  exp - s-  (  £  u.  G(tv+s))J 

V+  L  2  ks1  k  k  j 


4*oo 


=  exp  -  A  J  d|*(a)  /  ds  [l  -  exp  -  y  (  £  u1G(tv+s))al 

%  —  L  k=1  K  K  J 


because,  for  t^  +  3  <  0  or  s  <  -  and  so  for  s  <  0,  one  has 
G(tk  +  s)  *0. 

Furthermore  since  any  translation  on  s  does  not  change 
the  integral  the  assumption  tv  ^  ttv  tav  ...  v  t  *  0  is 
not  a  restriction.  ■ 


/ 
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Remark  :  The  process  Z  as  defined  is  real  valued.  It  is  however 
very  easy  to  define  a  vector  valued  process  as  follows  : 

+  0D 

(21)  Z(t)  *  l  C.  G ( t-t . )  6 

j  =  —  CD  J  ^ 

Where  G  is  a  function  with  values  in  &q. 

To  obtain  the  characteristic  function  of  Z,  it  is  useful 
to  introduce  some  new  notation. 

Let  u  be  the  (q*n)  matrix  with  entries  u.  ,  and  Z 
be  the  matrix  with  entries  Z^t^). 

Then,  <u,  Z  >  ^  trace  ufc  Z  is  the  Hilbert-Schmidt 

inner  product  for  matrices. 

One  can  then  state  : 

1 1 .  5  -  Tfieo>te/a  : 

If  Z  is  a  filtered  Poisson  process  with  values  in 

(22)  1  tfifc  Z(t)  =  l  C.  G(t  -  t.) 

3  ~  3 

and  if  /  ||G(s)  ||  ds  <  ®  and  J  a2  du(a)  <  «  then 

* 

*  m  *  \J(t! ,  t,  ,  . . . ,  t  )  t  \\ m  *  f  ff?xm 

m  ’  u  ^ 

a  i<u,  Z>  \ 

(23)  *(u)  S  e (e1  '  * 

exp  -  X  /  dsdu(a)  £l  -  exp-^-  (  £  <uk,G(tk+s) » *j 

Remarks  :  Formulae  23  and  14  yield  the  same  result  if  in  (14), 
uk  G^tic+8^  is  replaced  by 

<V  G(tk  +  s)>  *  .1,  a  ui,k  Gi(V9) 

with  the  convention  that  ||G(s)  ||  2  ^  <G,  G>  «  Trace  Gt  G. 


/ 
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III  -  CALCULATION  OF  THE  MOMENTS  OF  Z 


A  new  process  Z  with  characteristic  function 

♦  (»>  =  Z(V)  -  E(el<U'Z 


(D 


exp-X  f  ds  du (a)  M  -  exp  -  yy  (J  uk  G(tk+s))*] 


is  useful  if  it  is  different  from  the  usual  one s,  Poisson  proces¬ 
ses,  spherically  invariant  processes/ and  more  accurate  to  some 
applications. 

The  aim  of  this  section  is  to  obtain  the  moments  of  Z.  To  that 
end  it  is  assumed  that  Z  has  moments  of  all  orders  and  that  $ 
is  analytical  at  the  origin. 


.k 


t  r>  \ 


*  ("}  =  F I  l  /u.Z/t,)  ♦  ///  +  u  Z  (t  )  t) 

-  u  K !  *  *  '«  »• 


k=0 


in  ui 


=  l  E  (<U,  Zr>k) 


k=0 


% 


The  explicit  expression  for  41  shows  that  4  (u)  is  a  func¬ 
tion  of  the  square  of  <u,  G->(s)  defined  by  : 


(3)  <u ,  G_>  (s)  -  I  u.  G ( t.  +  s ) 


III. I  -  Lemma  :  The  odd  moments  of  Z  are  equals  to  zero. 


P*oo d  :  The  series  expansion  of  the  exponential  function  yields 


(4)  1  -  exp  -  <u ,  Gr>*  (s)  ■  —  l  (-  %-) 


ap  <u,  Gg>2p(s) 


p«1 


Pi 


so  that 


(5)  <Mu)  »  exp  ♦  X  /  dsdu(a)  \  (- 


a  p  <u /  Gg>2p(s) 


Pi 


p«1 


/ 
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Finally  this  last  exponential  can  be  expanded  to  yield 
a  new  series  in  the  variable  j  ds  <u,  Gg>*  (s)  where  each 

R  + 

term  is  homogeneous  in  u  and  of  even  degree. 

The  odd  moments  which  are  the  coefficients  accompanying  expres 

i 

sion  of  the  form  u}1  u*2  u  ra  with  ii  ♦  i2  +  •••  + 

in  *« 

vanish.  « 


It  is  now  convenient  to  introduce  the  second  characteristic 
function  p  (u)  : 

(6)  <|i  (u)  =  log  p  (u)  =  -  \f  dsdu(a)  1  -  exp  -  <u,  Gc>*  <SM 

>*  + 

1 1 1.2  -  Lemma  :  ^  (u)  *  -  £  — — ~ 

q=1  4 

pf looj  :  p(u)  =  log  <t>  (u)  =  log  [l  -  (1-*(u))]  and  for 
0  <  e  <  1 

r*  eq 

log  (1  -  e)  =  -  (e  +  -j*  +  •••  ♦  ■q"+***)  * 


111.3  -  Lemma  :  p (u) 


co  p  oo 

I  I 

1=1  q  Ht*1 


(-1) 


i+1 


(21)  ! 


E  ( <U ,  Z 


>21i 


1 


Pfioo_^  :  Since  the  odd  moments  of  Z^ vanish,  the  expansion  of  p 
yields  : 

(7)  *(u)  =  l  E  (<u,  Zf>2£)  so  that 

(8)  1  -  <Mu)  =  l  {-jYTn —  E  (<u»  zg>  >  ■ 

Comparing  the  two  expressions  for  i|»  given  respectively 
by  Lemma  III. 3  and  the  relation  (5),  one  obtains  the  moments 


of  Z.  One  has  : 


/ 


(9) 


-l  - 

L  n 


q*  i 


l 

i  =  '\ 


(-  D 

UT) ! 


4  +  1 


E(<u,Z„>2*) 


*  <l<(u)  = 


X  l  ^7  /  du(a)  (-  ^-)P  /  ds  <u,G  >2p(s) 
p=1  P‘  IR+  1  1R  4 
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We  consider  first  the  coefficients  of  ul1  u22  with  ii  +  i2  =  2 
to  compute  moments  of  second  order. 


111.4  -  P?io position  :  Vt  E(Z(t)2)  =  A  {  a2du(a)  /  G(s)2  ds 

R+  IR 

^ti,  t2  £  (R  E  (Z  (t  i )  Z  (t2 ) )  =  a/  a*du(a)  /  G  (t2  +s)  G(t2+s)  ds 

IR+  IR 


Pa.0  0  ^  :  One  must  have 

1+1 


1  (-1) 
1  2! 


E  (<u,  Zr> 2)  = 


X  X-  f  du(a)  (-  %-  )  f  ds  <u,  Gr>2  (s) 


Using  the  respective  definitions  of  <u,  Z^>  (II. 15)  and 
<u,  G^>  (3)  one  obtains 

(10)  EtUiZttJ  +  u2  Z  (tj ) ) 2  = 

X  J  a*du(a)  j  ds  (UjGtti+s)  +  uaG(ta+s))  and  equating 

nu  * 

tht;  coefficients  of  uj  ,  uj  and  UiUa  : 


E (Z (t i ) * )  =  x/  a*du(a)  J  ds  Gttx+s)2  * 

KU  * 

Aj  a2 du (a)  J  ds  G(s)2  *  E(Z(t2)2) 
and 

E (Z (t x )  Z(ta))  *  X  /  a2du(a)  /  ds  G(tx+s)  G(ta+s)  • 

ft 


Remark  : 


If  ti  <  t*  ,  the  integrals  /  ds  G(ti+s)  G( t, +s)  and 

CO  R 

/  ds  G(s)  G(t2-ti+s)  are  equal. 

0 


The  next  proposition  lists  the  moments  of  fourth  order. 


HI. 5  -  P^oppA-tt-cow 

1)  E(Z(t4))=3x  /  a-du(a)  /  dsG4 (s) +3X2 


IR, 


IR. 


1 

/  a2du(a)  /  dsG2 (s) 


2)  E(Z(t1)2Z(t2)*)=3x/  a“du (a)  JdsG2  (t2 +s) G2  (t2 +s) 


IR. 


IR 


+  X2  [  J  a2du(a)  /  ds  G(s)2]2 
L  «+  J 

+  2  X2  /  a2du(a)  /  ds  G(ti+s)  G(t2+s)*|J 
TR+  IR  -1 

3)  E(Z  (ti)  Z  (t2)Z(t3)  Z  (tH)) 

*  3X  /  a“du(a)f  ds  G(ti+s)  G(t2+s)  G(t3+s)  G(t„+s) 

«+  * 

+  \  [  /  a2dw(a)j  If  ds  G(ti+s)G(tj  +s)  /  ds  G(tk+s)  G(tt+s) 


/ 
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where  the  summation  is  over  all  permutations  (i,  j,  k,  1)  of 
(1,  2,  3,  4)  such  that  i<  j  and  k  <  1. 


Pa.0  0 ^  :  The  first  terms  in  the  expansions  of  formula  (9) 
are  restectively ,  for  : 

(11) 


1 

(-D1+1 

E (<U,Z^>* ) 

.  (-ir  ' 

A  | 

E(<u,Zr>4) 

1 

2  ! 

O 

n  • 

1 

[(-d1+1 

E (<U,Z  >*) 

Z 

,  (-D2*1 

E(<u,Z^>4) 

J 

2 

2! 

4! 

+ 


Equating  the  homogeneous  terms  of  degree  4  in  u, 

(12)  E  ( u  2  Z  ( t  x )  +u2Z(t2)  +u3Z  (t3)  +u4Z  (t4)  )  4  + 

( 1 3 )  -  j  [jr  E<uiZ (tj  +u2z ( 1 2 )  +u3z (t3)  +u4z (tj ) *1 


one  obtains  : 


2 


1  ^ 

(14)  x  JT  /  X  du(a)  /  ds 


UiG (t i +s) +u2G(t2+s) +u3G(tj+s) +u4G (t4  +s) 


This  yields,  if  t .  =  t,  i  *  1 ,  2 ,  3 ,  4 
TT  ■<*<«*>  ■  T7-i  [»<*<«*]’ 

♦  X  2  ^  /  a"  d 4(a)  /  ds  G4  (s) 

The  first  term  on  the  right  hand  side  is  given  by 
Proposition  III. 4 


In  the  case  of  ti  ■  t3  <  t*  »  t* 


I 
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JT  C24  E(Z(ti)J  Z  (t  j  ) 1 )  = 

\  (^)  ^2E(Z(ti)J)E(Z(ta)J)  +  21 (E(Z (tx) Z (t2) ) ) J J  + 

1  2  £  4 

X  21  c-  du  (a)  /  ds  G1  ( 1 1  +s)  G*  (t2  +s) 


4 ! 

where  C*  =  — : — .  Again  the  first  term  on  the  right  hand  side 
(2 ! )  * 

is  given  by  Proposition  III. 4. 

The  case  of  tj  <  t,  <  t3  <  t4  is  more  complicated.  For  any 
random  variables  X±  and  real  number  u±,  i  =  1,  2,  3,  4,  one 
has  that  : 


(15)  [e(uxXx  +  u2X2  +  u3X3  ♦  u4X4) * 


2 


22  22  22  22 

♦  u».X;. )  ♦ 


2E(uiU2XiX2  ♦  U1U3X1X3  ♦  UiUHXiX4 


+  U2U3X2X3  ♦  U2U4X2X4 


+  u3u4X3X4) 


=  l  U l  [E(xiM]i 

*  4  jj  uiuj  [E,xi  v]' 

♦  2  [  u*u*  |e(X*)  E(X*)] 


* 4  i,A<,  uiuj  ukut  [e<xixj'  E*w] 
and  Tl'j.  5  (k,t. 


2  |  ui  uku«  [E,xi>  E(xkx«'] 


k<t 


/ 
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The  pair  [(i  <  j)  ,  (k  <  1)^  is  not  orderwlso  that  there  are 
6  «  5  such  pairs.  The  same  is  true  for  the  pairs 
[(i)  ,  (k  <  1)]  . 

One  then  equates  the  terms  containing  u^  u^  u^  u^  with 
(i,  j,  k,  1)  permutation  of  (1,  2,  3,  4)  and  i  <  j  <  k  <  1; 
there  are  only  six  such  permutations  and 

JT  41  E(Z(tx)  Z(t2)  Z ( t j )  Z ( 1 3 ) )  = 

I  E  ( Z  (t . )  Z(t.))  E(Z(t.)  Z(t.))  + 

i< j ,  k<l  1  3  k  x 

(i,j,k,l)  a  permutation 

1  4 

X  2  /  T"  /  ds  41  G(ti+s)  G(ta+s)  G (t 3  +s)  G (ti,  +s) 

**  A 

For  the  summation  £  one  choose  i  <  j  among  1 ,  2 ,  3 ,  4 

i< j ,k<l 

and  (i,j)7Uk,l) 

and  k  <  1  are  the  two  remaining  integers.  • 

Remark  :  One  could  also  compute  other  moments  such  as 
E (Z (t i )  Z(t2)2  Z(t3))  and  E(Z(ti)  Z (t2 ) 3 )  but  these  are  of  a 


lesser  interest. 


i 
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IV  -  IPENTI FICATION  OF  THE  GREEN  FUNCTION  G 

Given  some  experimental  data,  one  may  assess  the  adequacy 
of  the  model  and  in  case  it  is  satisfactory,  to  estimate  the 
paratemers  which  are  the  Green  function  G,  the  parameter  X  of 
the  original  Poisson  process  and  the  mixing  law  u  for  the  ampli¬ 
tude  C.  In  this  section,  an  estimation  procedure  for  G  is  pre¬ 
sented.  Estimation  of  X  and  4  is  the  subject  of  the  next  one. 

Statistical  estimation  based  on  one  trajectory  of  the 
observed  process  requires  that  the  process  be  strongly  statio¬ 
nary  which  is  the  case  for  the  process  Z.  Its  characteristic 
function  (theorem  I 1.4)  is  indeed  invariant  under  translation 
of  the  time  variables. 


The  covariance  r  of  X(t)  and  of  X2  (t)  are  stationary 

and  given  oy  tne  relation  : 

(1)  r ( t )  =  E(Z (t) Z (t-x) ) =  X  /  a2du(a)  J  G(s)  G(s+x)  ds 


(2)  r  2  (x)  *  E(Z2  (t)  Z2  (t-x)  )  = 


3X  J  a-du(a)  J  ds  G2  (s)  G2(s+x) 


IR  , 


IR 


xj  a2du (a)  /ds  G2 (s) 


L  ®. 


|R 


+  2  X 3 


j" J  a2du(a)  /  ds  G(s)  G(s+x)j 
^R  IR  J 


In  (1),  x  /  a2du(a)  is  a  normalization  factor  and, 

assuming  (11.12)  one  has  r(0)  =  X  /  a2du(a). 

♦ 

IV. I  -  Pfropoa ition  :  Let  g(f)  ^  /  e+is^  g(s)  ds  =f^(G).  Then 

—00 

f  dp 

(3)  r (t)  »  r (0)  /  e~isf  | g (f ) | 1 

£  7f 


/ 
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g  is  the  Fourier  transform  )'  of  G,  hence  the  notation 
g  = J‘ (C) .  It  is  square  integrable  and  vanishes  for  s  <  0  ; 
one  defines  (cf  [3],  p.  30)  its  analytical  exten¬ 
sion  g(f  +  ib)  by  the  relation 

g(f  +  ib)  =  /  eis(f+ib)  q(sj  ds  b  >  0 

g  is  said  to  belong  to 
the  Hardy  space  H2+. 

0  f 

?\oo ^  :  Taking  the  inverse  Fourier  transform  3*  ^  of  g,  one 
gets  : 

,  <f  00 

(4)  G (s)  =  i'  *1  (g)  (s)  =  /  e"isf  g(f)  df 

In  (1)  the  integral  in  s,  is  in  fact,  a  convolution  in  s,  let 

(5)  G* (r)  4  G (-  r)  then 

(6)  /  G (s)  G ( s+t ) ds  =  /  G ( s-t ) G ( t ) dr  =  /  G'(t-s)G(s)  ds  =  G'*  G(t 

SI  *  *  s 

With  self-evident  notation,  we  know  that  : 


G’*  G 
s 

-f-' 

<<F 

G'  . 

1'  G) 

But , 

by  definition 

g  (f ) 

4? 

G(f) 

+  0D 

-  I. 

elsf  G(s)  ds. 

and 

g(f>  = 

fo- 

(f)  = 

/ 

—  ao 

eisf 

G  ( — s ) 

ds ,  consequently 

(7) 

r  ( r )  = 

r  (0) 

G’ 

s  G 

(t)  • 

= 

r  (0) 

27  / 

+  00 

e 

■isf 

|g(f) I1 

df  . 

This  result  shows  in  particular  that  the  gain  |g(f)| 
of  the  causal  filter  G  can  be  estimated  when  an  estimate  of  the 
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covariance  r  is  available. 


It  is  a  standard  procedure  in  automatic  control  to 
recover  the  response  G  of  a  causal  filter  from  its  gain  |g|. 

There  is  however  a  difficulty  which  can  be  illustrated 
with  the  example  of  a  rational  filter  G  whose  Fourier  transform 


is  given  by  the  relation  : 

(f  -  hi )  (f  -  b2)  ...  (f  -  b  ) 


(8) 


g (f )  * 


(f  -  ai)  (f  -  a2)  ...  (f  -  ap) 


p,  q  €  W 


Since  G  is  causal,  g  is  entire  in  the  upper  half -plane 
+  so  that  the  poles  a\ ,  at  ,  . . . ,  a  of  g  belong  to  the  lower 
half  plane  . 

+  oo 

Since  G(s)  =  J-  (  e-i^s  g(f)  df  is  assumed  real  and 

2  7T  ' 

—  oo 

then 


mm  -  i  r+” 

-  2n  ' 

—  OO 


_ __  .  +  «  .  _ _ _ 

A  Q  -  _L  f  rr  /  \  rlr 

g  —  2n  >  •  3  '  “ 


(9) 


g(f)  =  g(-f) 


(9)  shows  that  b  -b  are  simultaneously  zeroes  of  g 
and  a  -  a  simultaneously  poles  of  g. 

Let  b^  be  a  zero  of  g  and  define  g*  by  the  relation  : 
.  (f  -  bi)  ...  (f  -  b  »  f  -  b. 

9''f>  4  if  -  t-r-lpl 


f  -  b. 

l 


g'  is  a  frequency  response  such  that  |g'(f) 
f  -  b. 


| g (f ) | .  If  f  is 


real , 


f  -  b, 


has  modulus  1  and  is  called  an  interior  function 


(cf  [3] ,  p.  36) . 

f  -  b.  f  +  b 

Multiplying  (8)  by  -* - r—  *  -  ,  one  obtains  the  frequency 

r  "  Di  f  +  b± 

response  of  a  causal  real  filter. 

Thus  | g |  does  not  uniquely  determine  G  and  the  question 


arises  as  to  the  best  choice  of  G 


/ 
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IV.  2  '  VeAinition  :  A  causal  filter  G  is  said  to  have  a  minimal 


phase  if  —  is  the  frequency  response  of  a  causal  filter. 

For  a  filter  with  minimal  phase,  the  zeroes  b  of  g  have 
to  be  in  the  lower  half  plane  and  the  following  result  obtains 


([3]  p.  37)  . 


IV.  3  -  :  The  gain  |g|  determines  the  causal,  phase 


minimal  filter  G. 


In  fact,  the  frequency  response  g  of  this  phase  mini¬ 
mal  filter  G  has  an  analytical  extension  in  the  upper  half 


plane  given  by 


(10)  g  (f +  ib)  =  exp  4^-  / 


1_  (•  y  (f+ib)  +  1  lo 

lit  '  Y"(f  +  ib) 


Y2  +  1 


This  exDression  mav  be  used  to  aooroximate  a( f )  , 


f  e  ik  ;  for 


log  g (f +ib)  =  -  / 


b  r  log  |g(Y) Idy 


(Y-f)1  ♦  b2 


1  f  [il — l) -  +  — 1—  log  |  g  (  y  )  I  dY 

'■(Y~f  )  *  +b2  Y*  +) 


and  thus 


i)  k  /  lQg  ig(v) 

n  ;  (Y~f) 2  +b2 


dY  —  if—  log  1  g  (f ) 


([3]  p.  37) 


ii)  i  /  I" - *-1X1 _  +  _X_1  log  |g(Y)|  dy  „  Arg  g(f) 

L(Y-f)2+b2  y2+1J 

([3]  p.  38  and  5*,) 

When  b  small  this  last  integral  gives  an  approximation 
of  the  phase  Arg  g(f)  of  the  uniquely  determined,  causal, 
minimal  phase  filter  G. 


G  can  be  also  approximated  directly.  For  b  >  0,  define 


G  by  the  relation 


/ 
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+  00 

(11)  Gb(t)  =  |  e"ift  g (f +ib)  df 


Since  g(.  +  ib)  converges  in  L1  (&  ,  df)  to  g(.)  as  b  goes 

to  zero  ([3]  p.  30)  and  since  the  Fourier  transform  is  an  isometry 
in  L2  (*,<&  ,  df) 


(12)  Gb(t)  =  27  /  dfe 


-ift 


exp 


r  1  f  1  (f+ib)  +  1  log  |g(v)  I  ]  dv 
In  >  y~  (f +ib)  1+1  J  dY 


converges  in  L2  (ft,  ,  df) 


to  G(t),  as  b  goes  to  zero, 


(13)  /  L2ixl  1  d y  exists  for  functions  g  in  H^+. 

Y2  +  1 

The  desired  approximation  Gb  of  G  requires  that  the  Fourier 
transform  of  a  complex  function  be  obtained. 


In  conclusion,  the  correlation  function  r (t )  determines 
the  causal,  phase  minimal  filter  G  which  can  be  estimated. 

The  same  is  true  for  r2  (x)  as  defined  in  (2)  : 


r2  ( x )  =  3X  J  a4d|i(a)  /  ds  G2 (s)  G2 (s+x) 
(K+ 

+  r(0)2  +  2r  (T) 2 


Up  to  a  normalization  factor  3X  /  a4du(a)»  one  has  an 

ft* 

estimate  of  f  ds  G* (s)  G2 (s+x) .  This  function  of  x,  also  yields 

V 

the  gain,  denoted  |g2 (.) |  of  the  causal  filter  G2  and  one  may 
compute  from  this  gain  |g2(.)|,  the  response  G2  of  a  causal, 
phase  minimal  up  to  a  normalization  factor. 

If  the  two  estimations  of  G2 ,  the  first  based  on  the  mo¬ 
ments  r (x)  *  E(Z(t)  Z(t-x))  and  the  second  based  on  the  moments 
r2  (t)  3  E(Z2  (t)  Z2  (t-x)) ,  are  equal^nearly  equa^  one  can  safely 
accept  that  G  is  indeed  a  deterministic  function. 


/ 
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V  -  IDENTIFICATION  OF  A  AND  u 

Here  the  Green  function  G  is  supposed  known  and  the  aim  is 
to  estimate  the  unknown  parameter  A  of  the  indurlying  Poisson 
process  as  well  as  the  mixing  law  u  determining  the  law  of  C. 


Having  assumed  in  part  III  that  the  characteristic  func¬ 
tion  ^  of  Z  can  be  expressed  as  a  known  series  whose  coefficients 
are  the  moments  of  Z,  the  relation 


00  P 

iMu)  =  Log  Mu)  =  l  u2p  —  x  / 

P=1  2P  p!  *  + 

should  yield,  after  proper  truncation, 
terms  of  the  moments  of  Z . 


a2p  d|i(a)  f  G ( s)  2p  ds 
ft 

estimate  of  X  and  4  in 


The  problem  can  thus  be  stated  as  that  of  determining 
x  and  4  given  the  equalities 

(1)  X  /  a2p  du (a)  =  a  p  =  1,  2,  ...»  i 

*.  p 

where  is  known,  is  obtained  by  dividing  the  coefficient 

of  u2p  in  the  expansion  of  <|*  by  (-  1)p  2P  pi  /  G(s)2p  ds 

ft 

Such  a  problem  can  be  restated  as  a  classical  moment 
problem  ([l]).  Indeed  if  u  is  symmetric  probability  measure  on 
(ft  ,01*)  such  that  4 j  4  the  problem  at  hand  becomes  that 
of  finding  4  such  that  i  k  €  N  : 

/  a2k  d4  =  /  a2k  d4  =  a. 
ft  ft+  K 

J  a  d4  =  0 
ft 

The  procedure  is  as  follows.  Orthonormalyse  the  sequence 
of  functions  {ak  ;  k  6  ftl  u  (0)}  in  the  Hilbert  space 
L* ,  4)  in  order  to  obtain  a  sequence  of  orthonormal  polyno¬ 
mials  {P^  »  kfcWu  {0}}. 


I 
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We  recall  the  results  about  these  orthogonal  Polynomials 
as  shown  in  [l] ,  chap.  I  : 

Let  : 


(2) 


(3) 


valid 


s^  -  J  a*  du (a)  moment  of  order  k  of  u 


50  Si 

51  s2 


s. 


k+1 


|  denoted  the  determinant 


k+1 


For  the  family  {P^ 


s2k 

kt  w  o  { 0 } } ,  these  relations  are 


(4)  Pk(a>  ah  »  1 


Si  S2 


Si  •  •  • 

sk 

S  2  •  •  • 

sk+i 

h+1 

h+k 

3  •  •  • 

a 

>  +  bkpk+i 

(a)  i 

and 


(6)  b_1  =  0 


(7)  ak  =  /  aiPv<a>l  2 

K  tR  * 


(8)  bk 


^ ^  .•  If  (i  is  a  symmetric  probability  measure 

on  )  then 

V  k  MW  V  a  (r  )X 

P2k  (-  a)  *  p2k(a) 


P2k+1 


P2k+1 


(a) 


(-  a) 


/ 
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P*oo :  PQ(a)  =  1  and  one  proceed  by  induction. 
Since  a  and  P2k(a)  are  orthogonal  by  construction 


a2k  =  l 

a 

>P2k(a) 

1*  du 

(a) 

=  0,  so  that 

a  P2k(a 

)  - 

b2k-1 

P 

r2k-1 

(a) 

rg 

SX 

li 

P2k+1 

p2k.1<- 

a) 

1 

b2k 

[-  3 

P2k 

(-  a)  - 

b2k-1 

Similarly  < 

a2k+1 

/  a 

tP2k+1 (a)^ 

*  du  (a 

a  P2k+1 

(a) 

3  b2k 

P2k(a 

)  ♦ 

b2k  P2k+2(a) 

P2k+2(_ 

a) 

3  P2k  + 

2(a) 

■ 

V  .2  -  CofiottaKii 


P2k(a)  = 


^  D2kD2k-1 


P2k+1 (a) 


V°2k+1°2k 


1 

S2  •  .  • 

s2k 

S2 

S  4  •  •  • 

s2k+2 

•  #  • 

1 

a* 

a2k 

1 

s  2  •  •  • 

s2k 

S2 

S  4  •  •  • 

s2k+2 

m  •  • 

a 

•  •  •  •  •  • 

a*  ... 

•  •  • 

2k+1 

a 

2k-1 

*2k+1 


a  P2k(a) 


P*ood  :  One  integrates  a  P2k  (a)  with  respect  to  u  and  clearly 

2k-1 

2k_  - - - 


P91p(a)  is  orthogonal  to  1,  a,  . ..,  a  and  similarly 

.2k 


P2k*1  is  orthogonal  to  1,  a . a  .  a 

Since  only  a  finite  number  of  moments  Is  avalaible,  u 
can  be  determined  if  it  is  a  finite  convex  combination  of  point 


masses 


32 


Let  thus  ii  be  defined  by  the  relation  : 

(9)  II  =  Hi  ♦  Hj  +  ...  +  u#6  where 

oi  02  i 

(10)  Ui=U»+...+U2®1  Hi  >  o  ...  4*  >  0. 

The  equalities  X  /  a^  du  *  o  then  become  ; 

>*  + 

X(ui  Oi*  ♦  4*  0J 2  ♦  ...  ♦  44  ot’)  =  ai 

X  (Ui  01J,J  +  4»  Oj*'*  +  •••  +  4.  O02**)  =  Q] 

(11)  *  * 

n  «*  n 

X(UQ  oi1  '3i  +  4a  ctj,’*S'+  ...  ♦  **4  )  =  a7 1 

There  are  as  many  equations  as  there  are  unknowns. 


V.3  -  Remg->tfe-i  :  When  u  is  symmetric,  the  odd  moments  vanish. 
One  writes  hi  equalities  but  2*.  comes  from  cdd  moment?  f-rnm 
symmetry . 

/  4 o-i  % 

The  last  equality  for  the  odd  moment  is  J  a  dti(a)  =  0. 

The  last  relation  X  j  a*4  du(a)  =  does  not  appear  in  the 
classical  moment  problem. 

We  propose  to  compute  u  from  the  21  first  relations  and 
after  to  obtain  X  from  the  last  one  : 

x  /  a44  du(a)  =  a _ 

In  fact,  we  will  eliminate  the  unknowns  u^, 
i  =  1,  2,  ...,  i  of  the  system  (10)  (11)  to  obtain  X  directly 
from  the  a  1 

Once  X  is  known,  one  has  the  usual  moment  problem. 


I 
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V .  4  -  PnjQ£0_6lti£n  :  Let  u 
s2k  =  /  a2k  dn(a) 


S2k+2.0  S2k+2 


2) 


i 

l 

i=1 

k  an  integer.  Then 


l  u  6  ;  and 

i  o^ 


’2k+2i 


S2k+2 . 1  S2k+4 


s2k+2d  S2k+2i 


s 


2k*2l+2 


s. 


2k+4i 

This  determinant  has  l  +  1  rows. 


P-noo  ^  :  The  proof  is  written  for  the  case  of  l  »  2.  We  start 


Ui 

+  ^2 

- 

1 

M.  l 

a!  + 

U>2 

2 

a  2 

=  S2 

Ui 

U 

°i  + 

u2 

u 

°2 

=  S  4 

Ul 

6 

Oi  + 

^2 

6 

°2 

=  s* 

These  equations  are  linear  in  hi,  u2 ,  the  first  two  determine 
m,  >*2  and  the  other  two  must  be  consistent. 


1 

1 

i 

Ol2 

02  2 

S  2 

=  0  and 

4 

4 

So 

0  1 

02 

Ol2 

02  2 

®  2 

1 

1 

S2 

oi  4 

02  4 

So 

*  Ol  0  2  2 

Ol2 

N 

N 

D 

So 

Ol6 

026 

S6 

Ol4 

0  2  4 

S6 

1 

i 

i 

Let  V  ^ 

: 

oi2  V,  A 

o22  and  S_  ^ 
o 

S2  Si  ft 

Ol  4 

0  2  4 

S» 

s2 

So 

s« 


I 
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So  and  Si  both  belong  to  the  span  of  Vi  and  Va  and  any 
other  vector  in  this  space  is  dependent  on  So  and  Si.  This  is 
the  case  of 


S2  - 


s„ 

s6 

So 


since  the  relation 


yjBi9  +  v 2 o i  8  —  s$  yields 


Ol  4 

024 

S  4 

1 

1 

S  4 

01  6 

02  6 

S  6 

—  4  _  4 

=  0  1  02 

Ol2 

022 

S6 

01  8 

O  2  8 

Sa 

0 1  4 

02  4 

se 

=  0 


So  finally 


i 

S  2 

S  4 

S2 

Si* 

Ss 

0 

s  4 

S6 

s8 

Starting  from  the  equations 

!2k 


2k  .  2k  _  o 

UlOi  +  U2<l2  -S. 


2k+2  .  „  2k+2  _ 

U  1  O  i  +P2O2  -s 


2k+2 


2k+8 
y  101 


2k+8 


2k+8 

A  similar  result  is  obtained,  that  is 


2k 

s2k+2 

s2k+4 

2k+2 

S2k+4 

s2k+6 

SB 

0 

2k+4 

s2k+6 

s2k+8 

The  generalization  for  arbitrary  t  is  obvious.* 
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We  next  show  that  A  which  appears  in  (11)  is  the  solution 
of  a  linear  equation. 

Indeed  by  dividing  (11),  by  X,  yields 

Ul+U2  +  -.-+UJl  =  1 

Ul01Z  +  U  2  G  2  2  +  ...  ♦  U^O^2  =  5Li.  4  S2 


4 1  ^  4i 

UlOl  +  W  2  0  2  +  ...  +  V 1° l 


and  by  proposition  (V.4) ,  one  has 


1 

Ol 

A 

Z 

A 

012 

a*  + 1 

A 

A 

A 

h 

a*  +  1 

llL 

\ 

\ 

\ 

A 

ai 

•  •  • 

ai 

Cl  2 

•  t  » 

“n+l 

ai+ 1 

•  •  • 

a2i 

41 


=  0 


0 


4 


s 


41 


or 


which  can  be  written  as 
a  linear  equation  as 
follows  : 


a2  «£+1 

0  a i  •  • • 

13)  A 

•  •  •  •  •  • 

=  - 

oi  o2  ...  a4+j 

“i+1  °2e 

•  •  •  •  •  •  • 

ai  ai  +  1  •**  a2t 

Rema ik*  :  If  u  is  as  above,  there  are  relations  between  the 
a^s  which  must  be  satisfied  because  of  proposition  V.4 

For  example,  one  must  have  : 


(14) 


a  i  a  2 

ai*1  al*2 


al+1 

°24«-1 


0 


/ 
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If  these  relations  are  not  satisfied  where  calculations 
with  experimental  data  are  carried  out,  one  must  chose  a  larger 
value  of  l. 

The  gaussian  case  is  associated  with  the  value  l  =  1 
and  the  equality  aio3  -  a,J  =  JaJduJ  a6du  -  /  a‘*du  =  0. 

Thus  one  must  compute  moments  of  order  six. 

When  u  is  a  convex  combination  of  point  masses,  the 
Hilbert  space  L2(^+#&  +  »  u)  can  be  characterized,  indeed 

A 

l/.  5  -  PjL££o*JJUon  :  If  u  is  a  convex  mixture  of  l  dirac,  the 
Hilbert  space  LM&+,(R  +  ,  u)  has  dimension  21. 


By  corollary  V.2,  one  has  an  explicit  expression  for 


I  »  f  K  t  J-l  T  <-»  S  1  f-  A  ^  a 

~  w  -j —  —  ■ —  ■  •  a » 


by  proposition  V.4 .  with  respect  to  the  measure  u,  one  has 
-^lE>2«.^a*  I*  d»i<a)  *  0  and  P2fc  is  linearily  dependent  of  1,  a, 

2i-1 

•  *  •  p  a  • 

The  same  occurs  for  P2t+i  and  more  generally  Ik  >  0  for 
P2t*k‘  • 


I 
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Thus  P-,  is  a  linear  combination  of  1,  P  ,  . ..,  P_ 

^  Z  Jt "  1 

but  P2|l  is  useful  for  the  identification  : 

:  The  roots  of  P2i  are  precisly 

±0!  ±  02  ,  - -  ±0 , 


P*ooJ  :  The  proof  is  again  written  for  the  case  of  i  =  2  ; 
one  has  : 


But  Pi,  (o  i )  =  I 


Pi(oi)  = 


1  s2 

s4 

s2  S4 

S6 

1  a2 

a4 

u  1  ♦  U  2 

Ui  0X2 +y2a2  2 

yxai2+u2CT2 

‘  yiai4+y2o24 

1 

a  1  2 

y  2 

y  2  a  2  2 

y  2  0  2 

y  2  0  2  2 

y  2  0  2  4 

y  2  0  2 

1 

ai  2 

01  4 

‘♦a...  -  *♦ 


.  6  +W  2  O  2  6 


so  that 


The  general  case  is  similar.  ■ 

Once  the  polynomials  P0,  Px,  . ..,  P2fc-1  are  known  the  values 
of  uk  follow  k  =  1,  . ..,  i  since  (  [l]  p.  22). 

1  ..  _  1 
1  uk  -  Ttt - 


l  |Pi  <ok)|’ 

i=0  1  K 


An  other  solution  to  compute  uk  k  *  1,  i,  t 
would  be  to  consider  the  l  first  equations  of  system  (10), 
(11)  since  \  and  o i ,  ...»  are  yet  computed. 
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In  conclusion,  when  it  is  assumed  that  C  has  a  spherically 
invariant  law  rather  than  a  Gaussian  law  and  when  the  mixing 
term  u  is  a  convex  combination  of  point  masses 

a  5 

u  *  i  u.  6  ,  we  can  adjust  the  moments  of  C  with  the  expe- 

k=1  K  °k 

rimental  moments  of  C  till  order  4i-1,  and  using  the  extra 
parameter  \  of  the  Poisson  process  equalities  hold  till  order 
4i.  The  extra  moments  of  order  4fc+k  are  determined  by  the 
shape  of  u  (see  proposition  V.4)  but  there  is  no  more  moments 
to  check  discrepancies  between  the  chosen  model  and  the  experi¬ 


mental  data. 
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